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1 Introduction 



The spectra of the highly charged ions (HCI) are currently under intensive 
experimental investigations [1] - [2]. Due to the strong field of the nucleus 
acting on the electrons in HCI, the relevant theory should be fully relativistic 
and based on the principles of QED. Up to now major experimental and 
theoretical efforts were concentrated on the evaluation of the energy level shifts 
due to radiative effects (Lamb shift) or to the interelectron interaction, on the 
hyperfine splitting and on the bound electron g-factor. The recent status of 
the problem was described in a series of books and reviews [3], [4], [5], [6]. 

Transition probabilities were studied less intensively, though a considerable 
amount of experimental [7] - [8] and theoretical [9] - [10] investigations also 
exists. The full QED theory for the transition probabilities was considered 
in [3], [11] on the basis of the S-matrix approach and in [12] employing the 
two-time Green function approach (see also [6]). 

In this paper we apply the most general QED line profile approach to the 
derivation of expressions for transitions probabilities in HCI with radiative, 
interelectron interactions and weak interaction corrections. 

The problem of the natural line profile in atomic physics was considered first 
in terms of quantum mechanics by Weisskopf and Wigner [13]. In terms of 
modern QED it was first formulated for one-electron atoms by Low [14]. 
In [14] the appearance of the Lorentz profile in the resonance approxima- 
tion within the framework of QED was described and the nonresonant cor- 
rections were estimated. Later the line profile QED theory was modified also 
for few-electron ions [15] (see also [16], [3]) and applied to the theory of 
overlapping resonances in two-electron HCI [17], [18]. Another application 
was provided to the theory of nonresonant corrections in HCI and in neutral 
hydrogen [19], [20], [21]. " 

It was found in [22] that the line profile approach provides a convenient tool for 
calculating energy corrections. The most natural way to calculate the energy 
shift in QED is to calculate the resonance shift in some scattering process. 
This way clearly indicates the limits up to which the concept of the energy 
of an exited states has a physical meaning - that is the resonance approxi- 
mation. The exact theoretical value for the energy of the exited state defined, 
for example, by the Green function pole, can be compared directly with the 
measurable quantities only in the resonance approximation when the line pro- 
file is described by the two parameters: energy E and width T. Beyond this 
approximation the evaluation of E and T should be replaced by the evaluation 
of the line profile for the particular process. For applications of the line profile 
approach to the calculation of energy level corrections see [23], [24]. 
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In [15], [16] the line profile theory in QED was developed on the basis of 
the evolution operator 5(0, oo). Such a consideration, though most natural 
for describing spontaneous emission and absorption processes, faces serious 
difficulties due to unrenormalizability of the corresponding S-matrix elements 
< f\S(0,oo)\i >. The resonance approximation itself yields correct results 
[15], [16], but the introduction of radiative corrections leads to the appearance 
of unrenormalizable expressions. 

Therefore in this paper we will use only renormalizable "full" S-matrix ele- 
ments < f\S(oo, —oo)\i >. 

During the last decades the parity nonconservation (PNC) processes in HCI 
became a topic of permanent interest (up to now only from a theoretical point 
of view) [25] - [31], [32]. 

In neutral atoms, where successful experiments and accurate theoretical calcu- 
lations were performed, the main uncertainty is introduced by electron corre- 
lation effects arising due to the complicated atomic structure. In few-electron 
HCI electron correlation does not play any significant role what greatly sim- 
plifies the theoretical treatment. The major corrections to PNC matrix el- 
ements in HCI are the radiative ones. This allows, in principle, to test the 
Standard Model of electroweak interactions beyond the "tree" level. It is es- 
pecially important that the experiments with HCI would provide these tests 
in the presence of strong fields. 

Partly the electroweak radiative corrections in HCI were calculated in Ref. [32]. 
In this work we present the full treatment of the QED part of the electroweak 
corrections based on the line profile approach. According to [32] these correc- 
tions are dominant in HCI and, unlike the so called "oblique" corrections, are 
strongly field-dependent, i.e. different in HCI and neutral atoms. 

The paper is organized as follows. In section 2 we begin with the description 
of the process of resonant photon scattering on the "tree" level. In section 3 
electron self-energy insertions (SE) in the electron propagator within the reso- 
nance approximation are considered. It is shown that they improve the energy 
denominator of the photon emission (absorption) amplitude. The same is done 
for the vacuum polarization (VP) insertions. In section 4 the standard expres- 
sions for the Lorentz line profile in the case of the photon emission (absorption) 
process is derived; in section 5 it is shown in which way higher-order radiative 
corrections to the energy can be incorporated in the Lorentz denominator. 
Section 6 is devoted to the derivation of radiative (QED) corrections to the 
initial state in the photon emission amplitude. 

The derivation of QED corrections to the final state requires more refined 
considerations, that are contained in Section 7. Here the process of double 
resonance photon scattering is introduced. In sections 8, 9 the SE and VP 
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insertions in the different electron propagators in the Feynman diagram de- 
scribing the double photon scattering are considered. These insertions produce 
energy shifts to both resonant energy denominators of the two-photon emis- 
sion amplitude. In section 10 the Lorentz profile for the one-photon emission 
process in the case of an unstable finite state is derived from the two-photon 
emission amplitude. This derivation includes the integration over the frequency 
of one of the photons. The resulting expressions can also be used to include 
the Lamb shift of the final state in the one-photon Lorentz denominator. The 
width of the final ground state formally can be set equal to zero; this approach 
can be considered as a kind of regularization of the singular amplitudes. Sec- 
tion 11 is devoted to the derivation of QED corrections to the photon emission 
amplitude (final state). 

In section 12 the vertex correction to the photon emission amplitude is evalu- 
ated; it is shown, that the ultraviolet and infrared divergencies, contained in 
the vertex cancel with the corresponding divergencies in the so called " deriva- 
tive" corrections. 

QED corrections to the amplitude admixed by the effective PNC weak poten- 
tial are investigated in Sections 13-15. It is supposed that the mixed states 
with opposite parity are separated only by the Lamb shift and it is shown how 
the vanishing denominators in the approximation of noninteracting electrons 
acquire a nonvanishing value. The cancellation of the divergencies in the vertex 
and derivative graphs with the effective PNC potential is also demonstrated. 

In sections 16-19 the results obtained earlier are generalized to few-electron 
ions. The first-order interelectron interaction is taken into account explicitly. 
The introduction of higher-order interelectron interaction corrections is dis- 
cussed. First-order interelectron interaction corrections to the PNC amplitude 
in the two-electron HCI are derived. Section 20 presents a short summary of 
all the results obtained in preceeding sections. 



2 Resonant scattering of a photon on an atomic electron 

In this section we consider first the process of elastic photon scattering on 
an atomic electron in lowest-order QED perturbation theory. As usual for 
bound-electron QED we employ the Furry picture in which the electrons are 
described by the solution of the Dirac equation 

0-m + 7 e\/)^ = 0. (1) 

Here jf = "y^p^,p k = —iV(k = 1, 2, 3),po = E, where m and e are the electron 
mass and charge, respectively, E is the energy, 7 M are the Dirac matrices and V 
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is the potential of the nucleus (point-like or extended). The pseudoeuclidean 

metric (H ) in 4-space is used. We employ relativistic units H — c — 1 

throughout the paper. The Feynman graph corresponding to the process under 
consideration is presented in Fig. 1. We will consider the elastic scattering and 
will assume that the initial (final) state A is the ground state. According to 
the standard correspondence rules (see, for example, [3]) the S-matrix element 
for the graph Fig. 1 is: 

= e 2 j d A x l d A x 2 ^ A {x 1 ) 1 , 1 A^S{x l) x 2 ) 1 ^J{x 2 )^ A {x 2 ) , (2) 

where ^ a(x) = ^ A(x)e~ tEAt is the electron wave function, ^! a is the Dirac 
conjugated wave function, A^*{x) = A^*{x) e~ luJt is the wave function (elec- 
tromagnetic potential) for an absorbed photon with frequency to. The electron 
propagator for bound electrons is taken in the form: 



oo 



S(x u x 2 ) = -1- e **<*-*> V duJl , ( 3 ) 



-oo 



where the sum runs over the whole Dirac spectrum. Integrating over time and 
frequency variables ti,t 2 ,Ui and using the relation between the S-matrix and 
the amplitude U 

S i f = -2ni6(E i -E f )U if (4) 
we will obtain for the scattering amplitude 

n E ni — Ea~ CO 

with the condition to = to', which implies energy conservation. Here we abbre- 
viate 

U^x) = e^A^x) . (6) 

In the resonance approximation the photon frequency to is close to the energy 
difference of two atomic levels: to ~ Ea> — Ea- Accordingly, we have to retain 
only one term in the sum over n in Eq. (5) 

jj(2)tcs = {U^)aA'{Uuj')a'A ^ 

Ea' — Ea — uj 

Eq. (7) shows that in the resonance approximation the scattering amplitude 
factorizes into an emission and an absorption part. It follows from Eq. (6) that 
the emission amplitude can be expressed as 

jj cm = ( U ^ AA ' (g) 
E A > — Ea — to 
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For the absorption we may write 



U ah = (Ma ' A . (9) 
E A , -E A -u> K ' 

This implies that the Lorentz profiles for emission and for absorption will be 
identical as it should be. Eq. (7) also indicates the limit, up to which the 
absorption and emission processes are independent - that is the resonance 
approximation. Only within the resonance approximation the spontaneous 
emission rate is independent from the process how the exited state is prepared. 

This approximation is quite adequate in atomic physics. Nonresonant correc- 
tions to the transition frequencies are still too small to be observed in present 
modern experiments [19] - [21]. 



3 Radiative insertions in the electron propagator 

In this section we consider radiative insertions in the electron propagator. 
The electron self-energy insertion is depicted in Fig. 2. The S-matrix element 
corresponding to this diagram can be expressed as 

S 4 = e 4 J d 4 Xl . . .d 4 x 4 ^ A (x 1 )^ 1 A^(x 1 )S(x 1 ,x 2 ) 

x ^ 2 S(x 2 ,x 3 )^ 3 D^ 3 (x 2 ,x 3 )S(x 3 ,x 4 )^ 4 A^(x 4 )^f A (x 4 ) , (10) 

where Z) WM2 (xi, x 2 ) denotes the photon propagator in Feynman gauge 

OO 

D^ 2 (x u x 2 ) = — : - / ^gMti-WM™ (11) 

2rci rio J 



-oo 



with r 12 = \x\ — x 2 \. Integrating over time and frequency variables and using 
again Eq. (4), we obtain the following expression for the correction to the 
scattering amplitude 

ni,n 3 ( E m ~ E A ~ w){E nz - E A - u) 



where £(£") is the electron self-energy operator defined by its matrix elements 

? 2 /^/ V2 

Ml I 

n 2 v ^12 ' nin 2 n2n 3 



(HE)Un 3 = £, E C^Ie^-e^)) , (13) 

no \ '12 /tiitioiiim 
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together with 



In the resonance approximation we have ri\ = 113 = A' and the correction to 
the scattering amplitude takes the form 

rr(4)res _ rr 2 ires 

sc " sc -E A -a;)' 1 J 

Repeating the insertions in the resonance approximation (the next term of 
this series is shown in Fig. 3) leads to a geometric progression. Resummation 
of this progression yields [14] 

rrrcs 

sc {E A .-E A - U y [ } 

where 

E A , = E A + (t(E A + u)) A , A ,. (17) 



Accordingly, in the resonance approximation the emission amplitude is repre- 
sented by the expression 

~ Ea' — E a —u> (18) 

The operator T<(E A + uj) may be expanded into a Taylor series around the 
value E A + oj = E A > : 

±{E A + u) = ±{E A ,) + (E A i — E A — cu)J:'(E ai ) + ..., (19) 

where Ti'(E A /) = ( 9 ) e=e a , ■ The first two terms of the expansion (19) 
are ultraviolet divergent and have to be renormalized. The renormalization 
of the first term to all orders in aZ is well known. A commonly employed 
covariant procedure has been developed by Mohr [33] (see further develop- 
ments in [4]). Within this approach the divergencies are canceled analytically. 
Originally Mohr's method was first applicable numerically to high and inter- 
mediate Z values, but recently the high accurate extension for Z — 1 was 
also carried out [34]. Another noncovariant approach is based on the partial 
wave expansion of S together with a numerical cancellation of divergencies 
for each partial wave has been proposed in [35], [36]. This approach is known 
as partial wave renormalization (PWR). A variant of the PWR where the di- 
vergencies are canceled partly analytically and partly numerically utilizes the 
multicommutator expansion method [37], [38]. 
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Apart from self energy (SE) also the vacuum-polarization (VP) insertions 
in the electron propagator of Fig. 1 should be considered to all orders in 
resonance approximation. The lowest order VP insertion is depicted in Fig. 4a. 
For the VP correction the Uehling approximation is frequently applied. In 
this approximation the bound-electron propagator in the VP loop (Fig. 4a) is 
expanded in powers of the nuclear potential ( in powers of aZ) and only the 
first term of the expansion is retained (Fig. 4b). The Uehling approximation 
leads already to fairly good results even for high Z values [4] . We will adopt this 
approximation in our derivations for simplicity. Accordingly, the VP insertions 
leads to the following modification of Eq. (17): 

E A , = E A + {t{E A + cu) W + {Vu)a>a> , (20) 
where Vu is the Uehling potential [39]. 



4 Line profile for the emission process 



In order to obtain the line profile for the emission process we retain the first 
term of the Taylor expansion (19) and consider the energy denominator in 
Eq. (18) with a shifted energy value 

E A i = E A i + (T,(E A i)) A ' A ' + (Vu)a'A' , (21) 

where 

(E(Ea>))a'A' = L A ® — -F A i . (22) 

L A ^ denotes the lowest-order electron self-energy contribution to the Lamb 
shift and T A / denotes the lowest order radiative width. The other lowest-order 
correction to the Lamb shift is the vacuum polarization 

(Vu)** = LIT . (23) 

Since the corresponding energy shift is real the vacuum polarization does not 
contribute to the width. To this approximation the emission amplitude reads 

U = (Kh* (24: ) 

&A> + L A' — &A — U) — A' 

where L A > = L^+L^?. A method to incorporate the Lamb shift and the finite 
width of the state A in the energy denominator of Eq. (24) will be discussed 
below in section 10. 

As the next steps of the calculation we have to take the square modulus of 
U cm and to integrate over the emission directions z7 of the photon and to sum 
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over the polarizations e. Taking into account the definition 

E / dV\{U:) AA ,\ 2 = T AA ,, (25) 



^rcs . 



where u; ros is the resonant photon frequency, V AA i is the partial width of the 
level A', associated with the transition A' — * A, we obtain an expression for 
the transition probability dW AA i(u) for the emission process 



dW AA <(uj) = -^-'Y^ J du\U cm \ 2 u; 2 du 

e 



1 T AA idu 
= 27 (E A , + L A i — E A — oof + \T\, ■ (26) 

Eq. (26) defines the Lorentz profile for the emission spectral line. The values 
for the resonance frequency in zeroth-order and in first-order approximation 
read: 

< s = E A , - E A , (27) 
ul s = E A , + L A i - E A . (28) 

For further derivations it is more convenient to introduce a complex resonance 
frequency 

ul es = E A t + {t{E A ,)) A , A , - E A . (29) 



5 Higher-order radiative corrections to the energy 



In this section we explain how to incorporate higher-order corrections to the 
energy levels in the denominator of Eq. (26). First, we have to take into ac- 
count the next term of the Taylor expansion (19). Considering this term as 
a correction to the energy denominator, we have to insert the resonance fre- 
quency given by Eq. (29). This leads to the "reducible" part of the second-order 
loop-after- loop self-energy (SESE) correction [3], [40]. 

A£ A ,(SESE,lal,red) = (£(E A ,)) A , A ,(p(E A ,)) A , A , . (30) 

The irreducible part of the SESE loop-after-loop correction can be obtained 
according to Fig. 3, where n 2 ^ A'. The corresponding contribution to the 
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scattering amplitude is given by 



_ (U^)AA' (E(E A + Uj))A'n 2 (X'(E A + Ld)) n2 A> 

E A ' — E A — to n ^ A , E n2 — E A — to 

^ <U u )a>a. (31) 

In the sum over n 2 7^ A' we can set E A + uj = E A >. This sum can be viewed as 
a complicated higher-order insertion in the electron propagator in Fig. 1. Re- 
peating this insertion subsequently in the resonance approximation, we finally 
obtain the additional energy shift of the denominator of the emission ampli- 
tude in Eq. (18). This shift represents the irreducible loop-after-loop SESE 
correction to the energy of the level A': 

A£ A ,(SESE, lal, irr) = £ (^OW^aOW (32) 

U2^A' "2 A 

The corrections (31) and (32) represent only two of the second-order radiative 
corrections to the energy. The full set of these corrections, i.e. the remaining 
second-order self-energy (SESE) corrections, the second-order vacuum polar- 
ization (VPVP) and the mixed SEVP corrections can be included in the de- 
nominator of Eq. (18) and hence in Eq. (26) in the same way. The expressions 
for these corrections can be found in [4] and [40], respectively. 



6 Radiative corrections to the emission amplitude (initial state) 



SE corrections to the amplitude (initial state), i.e. the SE correction to the 
wave function A' in the expression for the amplitude follows according to the 
diagram in Fig. 2 for intermediate states n\ 7^ A'. In this case the correction 



to U sc takes the form 



= - E 

ni^A' 



(U*) Ani (£(E A + cu)) niA , (U u ) a >a 



Em — E A + to 



E A i — E A — lu 



(33) 



Here in the sum over n\ we can set E A + u = E A >. Comparing the expression 
(33) with Eqs. (7) and (8) we observe, that Eq. (33) represents the corrections 
to the matrix element (U*)aa' i n the emission amplitude. Assuming, that all 
the improvements (corrections to the energy of the level A') in the denominator 
are already performed, we now can rewrite formula (18) in the form 



jj \~ uj) AA' 

E A i — E A — to 



(34) 
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where 



K U Jaa> = 2^ p p • (35) 

The vacuum-polarization contribution can be introduced immediately into 
Eq. (35) 

/ rm V- (U^AnA&jEA'^mA' + {Vu) ni A'\ , , 

\ U oj)AA' — ^ p • 

m+A' E ni - h A , 

In order to obtain the correction to the state A in the matrix element (U*)aa' 
as well as the correction to the energy Ea in the denominator of Eq. (34), we 
have to consider the diagrams with radiative insertions in the upper external 
electron line. However, these graphs appear to be divergent. The singularity 
occurring can not be regularized in a direct manner. To avoid this technical 
difficulty we may consider the double-photon-resonance scattering process of 
an electron in the state Aq with resonance absorption into the state A and 
finally into the state A'. If A is the ground state, A is some fictitious state 
that plays the role of a regulator. However the final expression for the line 
profile for the emission transition probability A' — > A does not depend on A 
and contains only a dependence on the width Ta of the state A. This width 
can be set equal to zero at the end of the evaluations. Such a regularization 
program will be carried out in the following sections. 



7 Double resonant photon scattering 



According to the idea discussed at the end of the previous section, we now 
have to consider the process depicted in the Feynman diagram of Fig. 5. After 
the integration over time and frequency variables the scattering amplitude 
corresponding to the process depicted in Fig. 5 takes the form 

ni E ni - E Ao - C^O n 2 E n2 ~ E Aq ~ UJ - UJ 

sr^ (U lv ')n2n 3 (U j' o ) n3 Ao , s 

X X E nz ~ E Ao ~ ^0 - U + U? ' 

The energy conservation during this process is implemented by the condition 

uo + iq = uj' + uj'q (38) 
and the resonance frequencies are given by 
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u = ^o> = E A - E Ao , (39) 
u = J = E A i - E A . (40) 

In the resonance case we have to set n\ = A, n 2 = A', n 3 = A in Eq. (37), 
which yields 

sc ~~ (E A - E Ao - cu ) {E A , - E Ao -u -lu) (E A - E Ao - u' ) ' 1 ' 

In order to describe the line profile for double-photon emission we have to 
consider 

rj {UZo)aqa{UZ)aa' U9 s 

m (E A - E Ao - uj ) (E A > - E Ao - u - cj) ' ^ 

i.e. the amplitude for double photon emission in resonance approximation. 
An analogous expression can be derived for the amplitude of double photon 
absorption 

jj _ (Uw>)a'a{Uu> )aA , , 

ab ~ (E A , - E Ao -u' - u') (E A - E Ao -u' )- [ ' 



8 Radiative insertions in the central electron propagator 



We start with the graph depicted in Fig. 6. Having performed the integrations 
over time and frequency we find 

ni E ni - E Ao - u „ 2 „ 3 {E n2 - E Ao - u - uj){E n3 - E Ao - lo - to) 

x {Uu>')n3n4\Uu>' )n4A (4:4:) 

- E nA - E Ao -uj -u + uj r 

In the resonant case (ni = A,n 2 = n 3 = E A >,ii4 = A) the equations above 
yields 



£j(6)rcs _ (Uu )a a(Uu)aA' 

(E A — E Ao — uj )(E A t — E Ao —ujq — u) 

x (H E A U}))aA' (U^A'AjU^AAp ^ 

E A i — E Aq — LOo — lu E A — E Ao — lo' 

where the condition (38) has been employed in the last denominator. The 
resummation of the SE and VP insertions to all orders in perturbation theory 
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leads to the following expression for the emission amplitude 

U em = WJma{K)aa. (46) 

(E A - E Ao - lu )(E a ,(lj + u ) - E Ao -u -u) 

where 

E a ,(uj + lo ) = E A , + (±(E Ao + LUQ + u)) A >A> + (Vu)a>a> • (47) 

This result is similar to Eq. (20) which has been derived for the one-photon res- 
onance case. We can also repeat a similar derivation as for the result Eq. (34). 
For this purpose we should take n\ = A,ri2 ^ A',ns = A' ,ri4 = A in Fig. 6. 
Accordingly, we find for the emission amplitude 

u em = (u: )a o a(u:) aa , ^ 

(E A - E Ao - lj )(E a ,(lj + u ) - E Ao -u -u) 

where (U*) AA ; is defined by Eq. (36). Eq. (48), as well as Eq. (34) describes 
the radiative correction to the emission amplitude that improves the wave 
function of the initial state A'. 



9 Radiative insertions in the upper electron propagator 



In this section we turn to the radiative insertions into the upper electron 
propagator of Fig. 5. To give an example the lowest-order SE insertion is 
depicted in Fig. 7. Performing again integrations over time and frequency in 
the corresponding S-matrix element yields 

_ ^2 [OW^Aq +^))mn 2 

n in2 ( E m - E Aq - UJ )(E n2 - E Aq - Ld ) 

^ (E n , - E Ao -ujq-uj) 2 ^ {E ni -E Ao -u -u + u') ' 

The resonant case is determinated by the conditions: n\ = A,n 2 = A, n 3 = 
A', n 4 = A, i.e. 



C/(6) 



res 

sc 



(U* ) AoA (t(E Ao + u )) AA 



E A — E Ao — u E A — E Ao — loq 



x (K)AA> {U^A'AjU^AAo ^ 

E A > — E Aq — luq — lu E a — E Ao — lo'q 

We can assume, that all the resonant radiative insertions into the central 
electron propagator in Fig. 5 are already taken into account. Repeating the 
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radiative insertions in the upper electron line in resonance approximation and 
summing up the resulting geometrical progression finally yields 

u em = -„ (u: o u oA (u:) AA , 

(E A (u ) - E Ao - u)(E A ,(u +u)- E Ao -u -u) 

together with 

E A (u ) = E A + (±(E Ao + uj )) aa + (Vu)aa ■ (52) 

Eq. (51) represents the expression for the double-photon amplitude where both 
energies E A and E A i are improved due to radiative corrections. 



10 Line profile for the photon emission with energy corrections to 
the final state 



For this purpose it is sufficient to keep only the leading terms in the Taylor 
expansion of the denominators in Eq. (51): 

E a (lu ) =E a + it A (E A )) AA + (Vu) AA = E A + L A — l -Y A , (53) 



E A ,(cu + uj) = E A , + (Z A ,(E A/ )) A , A , + (Vu)a>a> = E A , + L A , - -T A , . (54) 

Integrating over both photon directions and summing over the polarizations 
(see section 4) an expression for the transition probability [41] of double- 
photon emission is obtained 



dw A ,^ Ao - (2?r)2 {Ea + La _ e ^ _ ujo _ 

x F AA , dudu ^ 

\E A > + L A i — E Ao — u — to — \T A i\ 2 

where T AoA , T AA > denote the partial widths as defined in section 4. In a next 
step we integrate over ui Q in the complex plane. The integration can be ex- 
tended along the real axis — oo < uj < oo since only the residues at the poles 
are contributing. We choose the contour in the upper half-plane. The poles 
are located at 

) = ^AAo + L a + l -T A , (56) 
uJq 2) = uj A ' Ao + L A , - to + ^T A , (57) 
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where u>aa = Ea — Ea and u)a>a = Ea> — Ea , respectively. As the result of 
the integration a new expression for transition probability dW aa' is derived, 
where the Lamb shift and the width of the state A are taken into account 



dW A A' — — FaoaTaa' 
2tt 



1 



X 



Ta \^a'a + — loaaq — La — f (Ta' + Fa) — w 



r- 7^\du. (58) 

Ia' \u>aa + L a - uja'A - L A < - \\J. A , + ^a) + oj\ 2 J 

After simple algebraic transformations Eq. (58) can be cast into the form 

, w 1 ^AoA^AA' (Ta + r^') duj 

dW A A> = 77-— FT T~ \2 T i/r ^2 ( 59 J 

2vr T A r A / (cj am - u) 2 + ^{T A ' + TaY 

where Cja'a — Ea> + La> — Ea — La ■ The latter reveals how the Lamb shift La 
of the final state A enters the Lorentz denominator in the expression for the 
emission probability. 

For simplicity, we may assume that Ta a — Ta an d L 'aa' — Ta ; - This implies 
that both states A' and A have only one decay channel: A' — > A and A — > A . 
Finally, we arrive at 

dW aa' = TTT- \2 , Ur i p a 2 • 60 

2tt (wa'a - ^) 2 + i( r A' + Ta) 2 

This formula does not contain any dependence on the state A . This state 
enters only indirectly through the definition of Fa- HA denotes the ground 
state, we can set T A = in Eq. (60) which yields 



TA'duj 
2tt (u> A >A-w) 2 + \r 2 A , 



d^V aa' = -- ^r-TT^r ■ (61) 



The expression above deviates from Eq. (26) only by the presence of the 
Lamb shift La in the denominator. The introduction of the state A Q acts as a 
regularization and Ta plays the role of the regularization parameter. In case 
of the ground state A this parameter is removed by setting it equal to zero at 
the end of the derivation. 



11 Radiative corrections to the emission amplitude (final state) 



Now we are in the position to determine the final-state radiative correction 
to the emission amplitude, i.e. to correct the final-state wave function. For 
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this purpose we return to Fig. 7 and consider the case: n\ = A, ri2 ^ A, n 3 = 
A', n 4 = A. We assume that the resummation of all radiative insertions in 
resonance approximation in Fig. 7 is already performed. Instead of Eq. (51) 
we can write 

u em = -, (u: )a a(u:)aa' ; (62) 

(E A (u> ) - E Ao - uj )(E a ,(uj E Ao -uq-u) 

where 

, Tr¥ . \{^{E A )) An2 + (Vu) An2 (C/*)„ 2 a' 

(K)aa> = E ~ p p • ( 63 ) 

n 2 ^A E n2 - E A 

Comparing Eq. (63) with Eqs. (34) and (36) one can see that it represents the 
corrections to the final-state wave function in the expression for the emission 
amplitude. 

Collecting now all the various corrections to the amplitudes and to energy 
denominators given by Eqs. (36), (63), (20) and (52), respectively, we can 
finally write the generic equation (61) into the form (for the ground state A): 

dW AA , = ± . (64) 

For convenience we reintroduced the notation T A / = T AA > and denoted by 
T^, the expression for the width T A > (i.e. for the transition rate A' — > A) 
where both wave functions of the initial (A') and final (A) states are corrected 
according to the formulas (36) and (63). These corrections to the width are 
additive, so that 

T AA , = T AA , + T Aa , + T aA , (65) 

where T Aa , and T AA , denote the corrections to T AA > which contain the im- 
proved amplitudes (63) and (36). 

In principle, the width T A r in the denominator of Eq. (64) also should be 
replaced by T AA ,. This means, that a rigorous evaluation of the radiative cor- 
rections to the transition probability should include the second-order radia- 
tive corrections in the denominator of Eq. (64). The imaginary parts of these 
second-order corrections would determine exactly the value of T AA , [3]. 



12 Vertex correction 



A third type of QED corrections to the emission amplitude, which we have to 
consider are vertex corrections. For the evaluation of these corrections we can 
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return at first to the one-photon resonance picture (Eq. (18)). The SE vertex 
correction to the emission amplitude corresponds to the Feynman graph shown 
in Fig 8. In the resonance approximation this correction reads 

, WW 

E A ,-E A -u' 

where is the SE vertex that is defined directly by the upper part of the 
diagram in Fig. 8. The renormalization of may be performed in momentum 
space [40] 



A( 1)rcn (/,i>) = A( 1 )(i>', i 6)- 7 ,A( 1 ), (67) 



where 



In Feynman gauge the counterterm 7 M AW is ultraviolet and infrared divergent. 
This divergency cancels, if we take into account another contribution to the 
SE vertex that follows from the graph Fig. 2 and which has not yet been 
considered. Without changing the results of Section 3 we can take the first term 
of the geometric progression with SE corrections and multiply the expression 
without radiative corrections (i.e., the right-hand side of Eq. (18)) by this 
term: 

d _ (£(E A + u;)) A , A , 

U em + U cm ~ ~ £ ^ p p — ■ l by J 

ti A i — L A — UJ tii A' — J^A — W 



Here we employ the notation U° m for the emission amplitude Eq. (18) without 
radiative corrections and U^ m stands for "derivative" correction that will follow 
from Eq. (69). 

Using the Taylor expansion (19) for T,(E A + uj) and substituting the resonant 
value (29) for the frequency uj = uj^H = Ea> + £U' — Ea into the energy 
denominator in Eq. (69), we obtain 

(^K)aA' 



U° cm + Ut m = - ~ " 1 + &(E A ,)) A , A , , (70) 

— &A — U) 



respectively 



jjo +U d +J jV 
em 1 em 1 em 



i A A' 



Ea< 



E a 
1 



UJ 



Eai — Ea — UJ 



{^K)aa> + VAi^)aa] .(71) 
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The renormalized expression for £' is given by [40] 

S' rCn (i>) = S'(i>)-E', (72) 

where 



(73) 



and E(^») denotes the lowest-order electron self-energy operator in momentum 
space. The counterterm E' contains ultraviolet and infrared divergencies. Due 
to the Ward identity 

A« = -£' (74) 

all the divergencies in Eq. (71) cancel. The total renormalized expression 
Eq. (71) seems to be asymmetric with respect to A and A'. A symmetric 
result can be obtained from the two-photon resonance scattering. 

Consider now Eq. (51). Repeating similar derivations as have been performed 
above for Eq. (18), we obtain 



U° +U d 



i (u*) AoA (u:) AA , 



2 (E A (u ) - E Ao - u)(E A ,(u +u)- E Ao -ujq-uj) 

A( E * ^ 



E A i - E Aq - W 0rcs - Urcs + S' A 

E^/ 

\ E A — E Ao - Chores + ^ A 

together with the lowest-order resonance energies 



+ (i, I =3=S ■ (™) 



CO, 



CO, 



(0) 
Ores 

(1) 
Ores 



= E A 



= E A 



E A , 

E Ao + E^ 



CO 



(i) 



-^A' — -^a + £a' 



(76) 
(77) 
(78) 



We omit here VP contributions in the denominators of Eq. (75). The values 
H A = (Ti{E A )) AA and H A ' = (£(-Ea'))a'A' are considered as being renormalized 



i (ren) 



E5f n) . Substitution of Eqs. (76), (77), (78) into Eq. (75) 



yields a symmetrized expression for the sum: 



U° + U d 

em 1 em 



(u*)a oA (u:) 



AA' 



(E a (uj ) - E Ao - <jj)(E a <(ujq + to) - E Aq - u 



00 



X 1 + 



V(ren) 



+ £' 



(ren) 
A 



(79) 
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The cancellation of the divergencies with those in U^ m still takes place. 

Now the expression (65) requires to correct the width by replacing V by 
V AA , where 

T AA , = T AA , + T Aa , + T aA , . (80) 

The width T AA > results from Eq. (26) by replacing the amplitudes U cm by the 
modified expressions Eq. (71) and (79), respectively. 

What further remains is to include the contribution of the VP vertex as de- 
picted in Fig. 9. In Uehling approximation the renormalization of the corre- 
sponding correction is straightforward. One has to replace the expression for 
the electron loop in Fig. 9 b), i.e. the photon self energy U(xix 2 ) by the known 
renormalized expression 11^ [41]. 



13 Weak interaction mixed amplitude 



To begin the investigation of additional corrections to the amplitude due to 
weak interaction we return to the generic one-photon scattering Feynman 
diagram Fig. 1 by adding at first the effective weak-interaction potential in the 
electron propagator (see Fig. 10). We can specify the corresponding scattering 
amplitude as 

jj = ^ (70^PNc)nin 2 (£4/)n 2 A 

m E ni — E A — lo n2 E n2 — E A — u 
together with the effective parity-nonconserving (PNC) potential 

G F 

V PNC = -^-^QwPn(t)75 ■ (82) 

Q w denotes the "weak charge" of the nucleus 

Q w = Z{1 - 4 sin 2 e w )-N, (83) 

which relates the numbers of protons Z and neutrons iV in the nucleus and 
the Weinberg angle Q w - the free parameter of the electroweak theory. The 
latter is usually determined via a comparison between theoretical predictions 
and data from various experiments both in high-energy and atomic physics. 
The recent adopted value for the Weinberg angle is <d w = 0.2312 [42]. 

Within the resonance approximation we set n 2 = A' in Eq. (81). Since the 
operator Vpnc has no diagonal matrix elements the term ri\ = A' is absent. 
However, we can choose a level ri\ = A" which gives a dominant contribution 
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to the sum over ri\. A standard example in atomic physics is: A = Is, A' = 
2s, A" = 2p. Accordingly in the sum over ri\ the term yielding the small 
denominator E A n — E& (in the resonance approximation) dominates. The 
emission amplitude reads 

tt ( U w)aa" (7oVpnc)a"A' (qa , 

Uem ~ p p p p • 1 84 J 

Hi A" — tii A — UJ tLi a' — Hi a — U) 

Introducing the resonance frequency in the first denominator we can further 
write 

jj = (Pu)aa» (7o^pnc)a"A' ^ 
E a" — E a 1 E a 1 — E a — UJ 

In particular for states A' = 2s, A" = 2p in the approximation of noninter- 
acting electrons Ea" = Ea>- This implies that in the first denominator in Eq. 
(85) one has to account for the Lamb shift AEl- In the next section we will 
show how it arises and how to use the form (84) for U em for this purpose. 



14 Radiative corrections to the weak interaction mixed amplitude 



The first step consist in the evaluation of the SE insertions in the upper elec- 
tron propagator in the diagram of Fig. 10 within the resonance approximation 
(see Fig. 11 for the first term of this sequence). 

The scattering amplitude, corresponding to the diagram of Fig. 11 is given by 



JJ =Y] (^)Anx yv (E(E A + Uj)) nin2 (-foVp K c)n 2 n 3 (Uuj')n 3 A ^ 

m Em - E a -uj „ 2 „ 3 (E n2 - E A - u)(E n3 - E a -uj) 

The resonant state is = Ea>- The states yielding dominant contributions 
are ri\ = n 2 = E A n. Keeping only these terms in the summations we obtain 

jj = ( U w)aa" (gCgU + u})) A "A" (iqVpnc)a"a> ^ g7 ^ 
Ea" — Ea — u} (Ea" — Ea — uS) Ea 1 — Ea — 

In resonance approximation we can set Ea + uj = Ea> in the matrix element 
(E(Ea + oj))a"A" and, furthermore, E A > = Ea>> due to the almost degeneracy 
of the levels A 1 and A". Then the resummation of the SE insertions yields 

T r _ (U*)aa" (7o^pnc)a"A' 



Ea" — Ea — Ea> — Ea — uj 

If rii = A" but n,2 7^ A", we obtain the correction to the PNC matrix element 
(Vpnc)a"a' (assuming that all insertions in the resonant approximation are 
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already performed) 

rr _ ( U u)aA" (10VPNC)A"A' /oqN 

Uem ~ ~~ " p p , l»yj 

ti a;i — hi a — OJ hi A! — Cja — IjJ 

VYoVpkc)a»a' = 1^ B B ' 

n 2 ^A" h n 2 - E A » 

together with Ea> = Ea" in the denominator. This represents the correction 
to the wave function A" in the PNC matrix element. 

In order to account for vacuum-polarization insertions we simply have to re- 
place Eq. (90) by 

\(E(EA"))A"n 2 + (Vu)A"n2 (7oVpNc)n 2 A' 

(7o^pncW = L B B — • ( 91 ) 



n 2 ^A» E «2 - ^A" 

A vacuum-polarization correction of this kind has been calculated numerically 
in [32]. 

Now we turn to the evaluation of the SE insertions in the lower electron 
propagator in Fig. 10 (see Fig. 12). The corresponding scattering amplitude is 



jj _ {Uw) Am (70^PNc)nin 2 (^(-^A + ^))n 2 n 3 (Uuj')n s A 

SC ~ ^E ni -E A -u^ {E n2 -E A -u){E n3 -E A -u) 

(92) 

In this case the dominant state is n\ = A" and the resonant states are = 
n 3 = A', respectively. The expression for the emission amplitude can now be 
expressed as 

jj _ (U^)aa" (7o^pnc)a"A' (£>(Ea + u)) A 'A' 
E a" — E a — u} E a' — Ea — uj \Ea' — Ea — uj) 

After resummation of SE insertions for the resonant state A' to all orders, we 
obtain 

jj _ (U^)aa" (7o^pnc)a"A' _ {Uu)aa" (iqVpnc) a"A' 
E a" — Ea — u> Ea 1 — Ea — uj E ' a" — Ea' Ea' — Ea — u> 

Here we have used the resonance condition uj = Ea> — Ea- 

Again, the summation of all SE corrections in Fig. 11 is supposed to be already 
performed. Moreover, we will assume that VP contributions are included in 
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both Fig. 11 and Fig. 12, respectively. Accordingly, we can write 

Um - E A ,-E A -u AE L ' (%) 

together with AE L = L A n — L A i. Eq. (95) holds for the degenerated case, when 
E A t = E A n. 

Setting ri2 ^ A' in Fig. 12 we obtain the correction to the PNC matrix element 
(respectively to the wave function of the state A"): 

U cm = - iU:)AA " (7oV 7 cW ' , (96) 

E A ,-E A -u AEl y ' 

where 

, Tr v-^ (70^PNc)A"n 2 (E(E A »))n 2 A> + {Vu)n 2 A> .. 

(To^pncW = E p ~ ( 97 ) 

A similar replacement i<M» — -em — ^ — > AE L can be performed in Eq. (89). 

Graphs with Vpnc inserted in the outer electron lines in Fig. 10 are less im- 
portant since they do not contain dominant terms with small denominators. 



15 Electromagnetic vertex correction to the weak interaction mixed 
amplitude 



The weak interaction amplitude with one-loop SE corrections at the PNC 
vertex is described by the Feynman diagram in Fig. 13. The corresponding 
expression for this correction, denoted as U^ c w , can be written as 

jjVW _ _ (^qj)Arei(Aj^ B IJ ') nin2 (U ul ,) n2A 

nin 2 ( E m E A w){E n2 — E A — to) 

Here A^ 1 ^ is the weak interaction vertex that differs from the pure electro- 
magnetic vertex as discussed in section 12 by changing the vector-coupling 
7^ that enters in the expression for Aj^ to the pseudo-vector matrix 757^. 

A weak interaction potential according to Eq. (82) may be introduced by 

B, = -5, -^=Q wPN (r). (99) 

Evaluating Eq. (99) within the resonance approximation implies n 2 = A' to- 
gether with rii = A" as the dominant state. 
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Again we assume that the resummation of radiative insertions in both electron 
propagators in Fig. 10 is already performed: 



U. 



vw 



Ea> — Ea — U) 



AEr 



(100) 



The vertex A^ 1 ^ is ultraviolet and infrared divergent. The renormalization 
scheme can be taken over from the pure electromagnetic vertex: 



A (l)W (r en) =A (l)^ (/? 0_ 7|iA (l>^ 



(101) 
(102) 



In view of Eq. (68) we can write 



7 „A'«»' = -Ave* J J^JL^,*^ = 7 5 7,A W . (103) 

Similarly, the cancellation of the divergency in U™ occurs through the deriva- 
tive corrections. Consider once more diagrams with SE insertions of the type 
as depicted in Fig. 11. Without changing the results of section 14 we can pro- 
ceed with the first term of the geometric progression in the same manner as 
it was done in section 12. This yields 



U ow + U 

^ em 1 w em 



dW 



I A A" 



(To^e 



PNC) A" A 



E A > -Ea-uj AE l 
Combining Eq. (104) with Eq. (100) we have 



- (l + E^) . 



(104) 



dW 



^ cm 1 w cm 1 w em 



AA" 



Ea> — Ei 



AEr 



X 



(7oV r PNc)A^ + (Ai 1) ^ (rcn) ^)A^ + sif n) (7o^p NC )^| • (105) 



Since Vpnc — 75-^0 we find that the divergent part of the second term in the 
square brackets in Eq. (105) is — A( 1 )(7 Vpnc) j 4"A', while the divergent part 
of the third term in square brackets is — £^'(toVpnc)a"A'- These divergent 
parts cancel due to the Ward identity. The total expression can be rewritten in 
symmetric form when performing the same manipulations with the first term 
of the geometric progression in Fig. 12. As a result, we find 



x 



0w 



U. 



vw 



dW 



) A A" 



1 



E A < — Ea — u> AE L 

'(ren) 



(7oVpnc)a"A' + (tf )W{ien) B») A „A> + ^T n) + S £? n) ) 



(7o^pnc)a"A' • 
(106) 
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In principle, there exists also the VP electromagnetic vertex correction de- 
picted in Fig. 14. However, as it was shown in [32], this correction turns out 
to vanish within Uehling approximation. 

Finally, we can represent all pure electromagnetic radiative corrections to the 
PNC matrix element in the form 

{joVpnc)a"A' = ('JoVpnc)a"A' + (7o^pnc)a»a' + (1oVpnc)a"A' 

+ (tf Mim) B") A „ A , + \ (Sj cn) + Sj, cn) ) { 1q V pnc )a»a> • (107) 

The VP parts of the corrections (7oVpnc)a"A' an d (7o^pnc)a"A' have been 
calculated earlier in [32]. 

Apart from corrections to the weak interaction matrix element there exist also 
corrections to the emission matrix element in the weak interaction emission 
amplitude. They either result by setting rii ^ A", but n 2 = A" in Eq. (92) or 
by inserting radiative corrections into the upper external electron line in the 
diagram of Fig. 12. In the latter case one may turn over again to the double 
photon resonance picture described in sections 7-11. The other corrections to 
the emission matrix element are given by the Feynman graphs in Fig. 15. 
These corrections can be treated in the same manner as vertex corrections in 
the case of a pure electromagnetic amplitude (section 12). 

Note, that there exist unseparable radiative effects which cannot be adjusted 
to the weak interaction matrix element or to the emission matrix element in the 
weak interaction emission amplitude. These corrections are depicted in Fig. 16. 
However, they do not generate small energy denominators of the type (AE^)" 1 
and thus are neglible compared to the corrections (107). Nevertheless, they 
have to be taken into account if for some reasons the small denominator is 
absent. 



16 The interelectron interaction insertions in the electron propa- 
gators 

In this section we begin to investigate the line profiles and transition rates in 
two-electron ions. The same formalism will be applicable to few-electron ions 
as well. 

For highly-charged ions the lowest-order radiative corrections and the lowest- 
order interelectron interaction corrections (IIC) are approximately of the same 
magnitude. While the first ones are of the order a = e 2 the second ones are of 
the order 1/Z which is the same for aZ ~ 1. Since the lowest-order radiative 
corrections and IIC are additive the radiative effects can be taken over from 
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case of the one-electron ions (sections 2-12). What remains is to include the 
lowest order IIC, i.e. one-photon exchange corrections. 



The first-, second- and partly the third-order IIC to the energy levels have 
been calculated within the framework of QED during the last decade for two- 
and three-electron ions [43]- [47], [24]. In particular, the line-profile approach 
was used for the evaluation of the reducible parts of the two- and three-photon 
exchange corrections in [24], 



However, the interelectron interaction corrections to the transition probabili- 
ties have not been studied so thoroughly. In this and the following sections of 
our paper we present a full QED approach to the evaluation of the IIC based 
on the line profile theory. We will mainly concentrate on the first-order IIC, 
although the generalization to higher orders is straightforward. 

Consider the elastic photon scattering on a two-electron ion with the first-order 
IIC taken into account. We assume that scattering occurs at the electron in 
the state A, while the electron in the state B plays the role of a spectator. 
This process is depicted in Fig. 17. For the description of the interelectron 
interaction we use the Coulomb gauge, distinguishing the exchange by the 
Coulomb and Breit (transverse) photons. The contribution of the graph in 
Fig. 17 to the scattering amplitude is given by 



jj _ V^iW^-"! V ' <^ /n\n;ri2-Q\^ u>' /n\A (108) 



(u*U ni ^ (V C ) n\A 

i — E A — L 

where Vc is the Coulomb interaction 



m Em — E A — u n2 E n , 2 — E A — u 



2 

(VcWab = / tti(£i)tt£,(£ 2 )— V A (x 1 )-* B (x 2 )d 3 x 1 d 3 x 2 . (109) 



,2 

We adopt the notation 

{F)a'B'-,ab = {F)a'B'ab — (F)b'A'Ab , (HO) 

where F denotes an arbitrary two-electron operator. The second term in 
Eq. (110) takes into account the "exchange" part of the interaction. 

As in earlier cases we employ the resonance condition E A < = E A + u and set 
n i = n 2 = A' leading to the correction to the emission amplitude in resonance 
approximation 

Hi A 1 — & A — ^ & A 1 — Ej a — W 

The next iteration with respect to the Coulomb interaction Vc follows from 
the graph in Fig. 18. In this case the resonance approximation implies that 
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ni, n-4 = A and n 2 , n 3 e A', _B. The inclusion of all possible combinations leads 
to the correction 

jr (U?u)aa> [{Vc)a'B;A'b} 2 ni9 x 

cm E A ,-E A -u{E A ,-E A -uy [iiZ) 

An all-order resummation of the iterations results in a geometric progression 
and thus in an energy shift in the denominator 

rr {u:)aa> , gx 

hi A i + l±£j a , b — h, A — uj 

Here AE C A \ B is the first-order Coulomb IIC: 

A£^ = (Vc) A , B;A , B . (114) 

Exactly the same derivation can be repeated for the Breit interaction leading 
to the expression (see Fig. 19): 

jj = {Uw)aa> ^ 115 ^ 

E A i + AE\, B E A lo 

where 

AE\, B = AE c A } B + AEf B , (116) 

^E A } B = (Vb)a'B;A'B- (117) 

The matrix elements of the Breit interaction Vb are defined as: 

(Vb Wb = -e 2 f^l) , (118) 

V Ti2 / A 'BA'B 

{V b )baia>b = e 2 ( a ^- l e^-E B \r 12 \ (n9) 

V r 12 / Bi'yl'B 



Here , aj 2 denote the Dirac matrices which act on one-electron wave func- 
tions depending on spatial variables X\, x 2 and i — 1,2,3 . 

Moreover, we can perform simultaneously the resummation of the radiative 
insertions in the electron propagator. According to the derivations presented 
in section 3, this yields the result 

U = ( u Z)aa> (12Q) 

em E A , + AE\, B -E A -u' 

The emission line profile can now be expressed, as described in section 4. 

Following the derivations of section 5 we account for the higher-order IIC as 
well. For example the second-order Coulomb-Coulomb IIC can be obtained 
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from Fig. 18 for 711,714 = A' but n2,n^A',B. The part of this graph con- 
taining two Coulomb interactions can now be considered as a complicated 
insertion into the electron line corresponding to the state A. Repeating these 
insertions and summing up the resulting geometric progression we shall obtain 
an additional shift in the denominator that will represent exactly the second- 
order Coulomb-Coulomb IIC. Similar manipulations can be performed for the 
second-order Coulomb-Breit and Breit-Breit IIC. In the two latter cases re- 
ducible contributions also arise due the dependence of the IIC corrections on 
the photon frequency u [24], [48], 



These reducible contributions, i.e., the so-called reference- state contributions 
or derivative contributions can be obtained in the same manner as the re- 
ducible contributions discussed in section 5. 

The derivation of the IIC correction to the emission amplitude follows to the 
derivation of the radiative correction presented in section 6. Assuming now 
that in Fig. 17 and Fig. 19 ri\ 7^ A' but ri2 = A' and that all previous insertions 
are also performed as well, we obtain the generalization of Eq. (120) 

U cm = ^ ( U Z)aa> ( 121 ) 

Ea> + ^E\, B — Ea — oj 

Here 

(tt*\ - sr (^Wi [(VcWijBA' + (VbWi^a'] noo x 

m+A> U ni - E A > 

is the first-order IIC to the emission matrix element (correction to the wave 
function of the state A'). 

To obtain the shift of the energy Ea in the denominator of Eq. (121), as well 
as the correction to the wave function ^a in the emission matrix element we 
have to turn over again to the double resonance photon scattering picture as 
described in section 12. The following section will be devoted to this issue. 



17 Double photon scattering on two-electron ions 



The double photon resonant scattering process described in section 7 can 
be applied without any modifications to two-electron ions. The first-order 
IIC insertions into the central electron propagator in Fig. 5 will yield as was 
shown in section 8 for the radiative insertions the same result that was already 
obtained with the single photon scattering: the shift of the energy Ea> in 
the Lorentz denominator and the correction to the wave function \1/ a 1 in the 
emission matrix element. Instead of repeating this derivation we start with the 
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IIC insertions into the upper electron propagator in Fig. 5. These insertions to 
first order are shown in Fig. 20. Applying similar steps as described in section 
9, we derive an expression for the double photon emission amplitude 

jj ^ (Uw )a a{UZ,)aai , 123 s 

cm (E A + AE\ B - E Ao - u) (E A , + AE\, B - E Aq - oj - u) ' 

For obtaining Eq. (123) we have set n\ = n 2 = A,n 3 = A',n 4 = A in the 
diagrams Fig. 20 a)-d). We can also generalize Eqs. (60), (61) for the Lorentz 
line profile both for the excited and for the ground state A. The only difference 
to Eqs. (60), (61) is that now the frequency uj a > a also accounts for the IIC 

&A'A{B) = E A + AE AB — E A i — AE A r B . (124) 

Such interelectronic interaction corrections can be of any order 

oo 

AE ab = J2 AE ab (125) 

i=\ 

and the same holds true for AE A > B . 



Accordingly we can replace lv A ' A ( B ) by its experimental value ijJ e A F A r B y This 
is especially important for two-electron ions, where the energy differences 
between the levels within a multiplet can be deduced experimentally more 
accurately, than from theory. 

Finally, we can write down the correction to the emission matrix element (to 
the wave function of the state A), following derivations presented in section 
11 

u em = ( u : )[ u :)aa> (126) 

(E A + AE AB - E Ao - cu ) (E A > + AE A , B - E Ao - lu - uo) 

where 

/rr*x \-( V c)BA;Bn 2 + (^3) BA;Bn 2 ] (^)n 2 A' no7 \ 

\ u u))aa' — p p • y lzl ) 

n 2 ^A ^"2 - ^A 

For deriving Eq. (127) we have set n\ = A,n 3 = A' \n± = A, but now with 
n 2 ^ A in Figs. 20 a)-d). 

Again we can write down all the corrections to the transition probability 
A'B — > AB (or the partial width T AB . A / B ) in the form of Eq. (65). 

The first-order IIC to the transition probability can be always represented by 
the sum 



AB,A'B 



^aa> + T AA > + F AA , (128) 
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where ^f^ A , and T^,^ represent one-electron partial widths with improved 
amplitudes given by Eqs. (122) and (127), respectively. 



18 The interelectron interaction corrections to the weak interac- 
tion amplitude 



To evaluate the IIC corrections to the weak interaction amplitude we can 
take the sections 13-14 as a guide-line. The generalization of the diagrams in 
Fig. 10 to two-electron ions with IIC insertions into the electron propagators 
is depicted in Figs. 21 and 22. 

Consider first the IIC in Fig. 21. Here we take n 3 = Ea> (resonant case) and 
n\ — ni — Ea" (dominant state). Then the summation of IIC insertions to all 
orders will result in (compare with Eq. (88)) 

jj = (UZ)aa" (7o^pnc)a"^ ^ 12g ^ 

Ea" + AEa" b — Ea — oj Ea> — Ea — oj 

We assume that the radiative insertions are already performed. If n 2 ^ A" we 
obtain the correction to the PNC matrix element (to the wave function ^ a") 

U cm = , (7qVpncW (130) 

Ea" + AEa" b — Ea — oj E a 1 — E a — 

where 

/ T7 \ [(Vc)BA";Bn2 + {Vb) BA" ;Bn 2 ] (7o^PNc)n 2J 4' / 1Q1 x 

n 2 ^A" h n 2 ~ E A >> 

For the IIC in Fig. 22 we have to choose the dominant state n\ = A" and 
the resonant states n\ — n 2 — A'. Then the summation of IIC insertions to 
all orders yields (assuming again that the summation in Fig. 21 is already 
performed) 

jj = ( u Z)aa" (7o^pnc)a"A' ^ 132 ^ 

Ea" + AEa>>b ~ Ea — to Ea 1 + AEa'b ~ Ea — w 

As in section 14 we employ the resonance condition uo = E^ + AEa> b — Ea 
obtained from the second denominator in Eq. (132) and insert it into the first 
denominator: 

jj = { U u)aa» (7oVpnc)a"A' ^ 133 ^ 

Ea" + AEa"b ~ Ea> — AEa'b Ea 1 + AEa'b — Ea — oj 

The meaning of "dominant" state for two-electron ions needs to be clari- 
fied. These states should now be specified in a different manner compared to 
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one-electron ions. The most advantageous situation for the observation of the 
effects in two-electron ions occurs when two levels belonging to the same multi- 
plet but possessing opposite parity appear to be nearly degenerated. Unlike in 
the case of one-electron ions in two-electron ions this can happen only acciden- 
tally. The standard example is the crossing of the ls2s 1 S and ls2p 3 P levels 
that occurs for Z = 6 (carbon), Z = 64 (gadolinium) and Z = 92 (uranium) 
He-like ions [25]- [29]), [31]. 

Introducing the frequency 

&a"A'(b) = Ea" + AEa"b — Ea> — AEa'b (134) 

we may replace it then by the experimental value <^a^a'(b)- The exact theoreti- 
cal value for the splitting lja"A'(b) between the nearly degenerated levels of the 
two-electron ions can be also used. Now the energy difference Qa"A'(b) plays 
the role of the small denominator AEl and we can write down the analogue 
of Eq. (95) in the form 

jj = _ (Uu)aa" (jqVpnc) a"A' ^ 135 ^ 
Ea'(b) - E a -uj AE l 

where AE L = Cja"A'(b) and 

Ea'(b) = E A , + AEa'b ■ (136) 

Setting now n-i ^ A' in Fig. 22 we obtain the correction to the PNC matrix 
element (to the wave function in the state A') 

jj _ _ ( U u)aa» (7qVpnc)a"A' ( 137 ) 
Ea'(b) - E a -uj AE l 

where 

/ T/ n \T- (loVpNc) A"n 2 [(Vc)Bn 2 ;BA> + (V B ) Bn 2 ;BA>] , 1QQ x 

VYovpnc)a»A' = 1^ p B • ( 138 J 

n 2 ^A> E n2 - E A > 

The similar replacement E A "(b) ~ Ea — oo — > AE L can be made in Eq. (130). 

Again we do not consider the graphs with Vpnc inserted in the outer electron 
lines, since they do not contain any small denominators AE L . 

Finally, we end up with an expression for the PNC matrix element with all 
the IIC taken into account 

(1oVpnc)a>'A' = (loVpNc) A"A' + (7o^ / pnc) a" A' + (Jo^pnc) A" A' ' (139) 

where the corrected matrix elements (7oVpnc)a"A' an< l (7oVpnc)a"A' are de- 
fined by Eqs. (138) and (131), respectively. 
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19 The interelectron interaction corrections in the Dirac-Hartree- 
Fock approximation 



The QED theory of highly- charged ions, based on the Dirac-Hartree-Fock 
(DHF) taken as zero-order approximation has been considered in [3]. Due to 
the energy dependence in the Breit-interaction matrix elements (see Eq. (119)), 
the Breit-interaction cannot be included exactly in the zero-order DHF scheme. 
Therefore, we will consider here the DHF approximation restricting to the 
Coulomb interaction only and leaving the Breit interaction to be accounted 
perturbatively. Nevertheless, we may emphazise that the DHF method to- 
gether with an approximate inclusion of the Breit interaction is widely used 
and leads to very accurate results even for high-Z ions (see, for example [11]). 

In DHF approximation the one-electron equation (1) should be replaced by 

(ft-m + -f eV + -f eV UF )^ = (140) 

where the operator Vhf is defined as 

N 

Vk F (0/(0= E / dV*+(f ')V c (r,f')^ m (r')f(r) 

m=l J 

N 

- J2 / d 3 r^i(f')V c (f,f')f(n^m(r)- (141) 

m=l J 

In Eq. (141) ^ m (f ) are the solutions of the DHF equation with eigenvalues 
E m , f(r) is an arbitrary 4-component spinor function and the summation is 
extended to all occupied states of the iV-electron atom or ion. 

For highly- charged ions we can consider the term 7oVhf i n Eq. (140) as a 
perturbation. Then the first-order correction to the Coulomb wave function 
^/a (the solution of Eq. (1)) will read 

*\(? ) = E P^^f ) • (142) 



Using the definition of the operator Vhf, we obtain 

E n — Ea 



*\(r) = E E {Vc) ; m ^ n{f > . (143) 

m=l n^A 



In case of two-electron ions the sum over m in Eq. (143) reduces to one term 
m = B. Then, we obtain the corrections (U*) AA - t and (U*) AA , to the emission 
matrix element (Eqs. (122) and (127) only with the Coulomb potential Vc) if 
we replace the corrected wave functions ^ A and \I/ A , by the expression (143). 
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This means that when evaluating the emission matrix elements with DHF 
functions we will automatically take into account the corrections Eq. (122) 
and Eq. (127) with the Coulomb operator. The DHF functions include also 
higher-order Coulomb corrections which are negligible for high-Z values. The 
same holds true for the corrections to the PNC matrix elements as given by 
Eqs. (131) and (138), respectively. 



20 Conclusions 



Summarizing the results obtained in the present paper we can state that a 
most general QED approach has been elaborated for the evaluation of various 
spectroscopical properties of highly-charged ions: energy level shifts, transition 
rates and line profiles. The physical process of photon scattering on atomic 
electrons provides a proper basis for rigorous QED evaluations of all quantities 
being of current interest in studies of highly charged ions. 

The resonance approximation has been proven to be a powerful tool for ob- 
taining nonperturbative results such as the energy shifts in the Lorentz de- 
nominators. This is a most physical way to determine the atomic energy shifts. 

Rigorous QED expressions for all generic radiative corrections to the photon 
emission process are presented including simultaneously the radiative shift of 
the Lorentz denominator. Thus a full QED treatment of the Lorentz profile is 
achieved. The weak interaction-induced parity violating amplitudes are also 
treated within this approach and a closed expression for radiative corrections 
to the PNC matrix element has been derived. This allows for future accurate 
calculations of the PNC effect in highly charged ions. Furthermore, the theory 
is extended also to few-electron ions. 

QED and PNC corrections to the emission amplitude are presented for the 
most interesting case of single-excited states in two-electron ions. The renor- 
malization procedure for both pure QED radiative corrections and QED cor- 
rections in the presence of a parity-violating effective potential is discussed 
and the cancellation of the infrared divergencies is verified. 

The general QED approach and the formulas for the various QED corrections 
derived in the paper provide a firm theoretical basis for numerical calculations 
in the near future. 
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A 
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Fig. 1. Feynman graph describing the elastic photon scattering on an atomic elec- 
tron. The double solid line denotes the bound electron in the field of the nucleus. 
The wavy lines with the arrows denote the photon emission and absorption. The 
indices A, n correspond to the initial (final) and intermediate states. 



36 
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Fig. 2. Feynman graph corresponding to the electron self-energy insertions into the 
electron propagator in graph Fig. 1. The notations are the same as in Fig. 1. 
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Fig. 3. Feynman graph corresponding to double self-energy insertions in the reso- 
nance approximation. 
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Fig. 4. Vacuum polarization insertions in the electron propagator. Part a) corre- 
sponds to the exact expression and part b) describes the Uehling approximation. 
The ordinary solid line in Fig. 4b denotes the free electron propagation, the dashed 
line with the cross at the end denotes the interaction with the nuclear potential. 
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5. Feynman graph describing the elastic scattering of two photons on an atomic 
ron in the state Aq. 
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Fig. 6. Feynman graph corresponding to the electron self-energy insertion in the 
central electron propagator in Fig. 5. 
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Fig. 7. Feynman graph corresponding to self-energy insertions in the upper electron 
propagator in Fig. 5. 
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Fig. 8. The self-energy vertex correction to the transition amplitude A' — > A in the 
resonance approximation. 
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a) b) 

Fig. 9. The vacuum polarization vertex correction to the transition amplitude 
A' — > A in the resonance approximation: a) exact expression, b) Uehling approxi- 
mation. 



44 



A 



'X/X/X/lt/V 

ni 




A 



Fig. 10. Feynman graph that describes the elastic photon scattering on an atomic 
electron in the presence of weak interaction. The effective potential of weak inter- 
action between the electron and the nucleus is denoted by a "sea-gull". 
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Fig. 11. The self-energy insertion in the upper electron propagator in Fig. 10. 
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Fig. 12. The self-energy insertions in the lower electron propagator in Fig. 10 
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Fig. 13. Feynman graph describing the self-energy electromagnetic vertex correction 
to the weak amplitude in Fig. 10. 
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Fig. 14. Feynman graphs describing the vacuum polarization electromagnetic ver- 
tex correction to the weak amplitude in Fig. 10: a) exact expression b) Uehling 
approximation . 
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Fig. 15. Vertex corrections to the emission matrix element in the weak emission 
amplitude a) self-energy vertex correction b) vacuum polarization vertex correction 
in Uehling approximation. 
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Fig. 16. Unseparable radiative corrections to the weak emission amplitude. Figs. 16 
a), b) correspond to self-energy corrections, Figs. 16 c)-e) correspond to vacuum 
polarization corrections in Uehling approximation. 
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Fig. 17. The Feynman graphs that describe the elastic photon scattering on a 
two-electron ion in resonance approximation with first-order Coulomb interelectron 
interaction corrections taken into account. Here A, B denote one-electron states in 
the two-electron ion, A' is the resonant state. The dashed line corresponds to the 
Coulomb photon. The graph Fig. 17a) describes the "direct" term in the amplitude 
and the graph Fig. 17b) corresponds to the "exchange" term. The contributions 
of the "exchange" graphs always enter with an additional minus sign which traces 
back to the antisymmetrization of the two-electron wave function. 
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Fig. 18. The Feynman graphs that describe the exchange by two Coulomb photons 
between atomic electrons in the process of elastic photon scattering on an atom. 



53 



B 



B 



ni = A' 
n 2 = A 1 



CO 'W^X/V 



m = .4' 



n 2 = A' 



B A 

a) 



B A 

b) 



Fig. 19. Feynman graphs that describe the elastic photon scattering on the 
two-electron ion in resonance approximation with the first-order Breit interelec- 
tron interaction corrections taken into account. The wavy internal line denotes the 
Breit (transverse) photon. The other notations are the same as in Fig. 17. 
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Fig. 20. Feynman graphs that describe interelectron interaction first-order inser- 
tions: a), b) correspond to the Coulomb "direct" and "exchange" graphs; c), d) 
correspond to the Breit "direct" and "exchange" graphs. 
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Fig. 21. Feynman graphs that describe the weak interaction-admixed amplitude 
of photon scattering on two-electron ions. The graphs a), b) correspond to the 
interelectron interaction first-order Coulomb "direct" and "exchange" insertions, 
the graphs c), d) correspond to "direct" and "exchange" Breit insertions. 
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Fig. 22. Feynman graphs that describe the weak interaction-admixed amplitude 
of the photon scattering on two-electron ions. Unlike the graphs in Fig. 21 that 
correspond to interelectron interaction corrections insertions to the upper electron 
propagator in Fig. 10, Fig. 22 corresponds to the interelectron interaction insertions 
in the lower propagator. 
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